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Abstract. We investigate the global dynamics of the universe within the framework of the Interacting Dark Matter (IDM) 
scenario. Considering that the dark matter obeys the collisional Boltzmann equation, we can obtain analytical solutions of 
the global density evolution, which can accommodate an accelerated expansion, equivalent to either the quintessence or the 
standard A models. This is possible if there is a disequilibrium between the DM particle creation and annihilation processes 
with the former process dominating, which creates an effective source term with negative pressure. Comparing the predicted 
Hubble expansion of one of the IDM models (the simplest) with observational data, we find that the effective annihilation 
term is quite small, as suggested by various experiments. 
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INTRODUCTION 

The detailed analysis of the available high quality cosmological observations (JUSIHQl and references therein) have 
converged during the last decade towards a cosmic expansion history that involves a spatial flat geometry and a recent 
accelerating expansion of the universe. This expansion has been attributed to an energy component (the so called 
dark energy) with negative pressure which dominates the universe at late times and causes the observed accelerating 
expansion. The nature of the dark energy is still a mystery and it is one of the most fundamental current problems 
in physics and cosmology. Indeed, due to the absence of a physically well-established fundamental theory, there 
have been many theoretical speculations regardingthe nature of the above exotic dark energy (DE) among which a 
cosmological constant, scalar or vector fields (see Q5L 0, 0, [it |H [IH and references therein). 

Most of the recent papers in this kind of studies are based on the assumption that the DE evolves independently of the 
dark matter (DM). The unknown nature of both DM and DE implies that we can not preclude future surprises regarding 
the interactions in the dark sector. This is very important because interactions between the DM and quintessence could 
provide possible solutions to the cosmological coincidence problem. Recently, several papers have been published 
in this area |[l2[ [l3ll proposing that the DE and DM could be coupled, assuming also that there is only one type of 
non-interacting DM. 

However, there are other possibilities, (a) It is plausible that the dark matter is self-interacting (IDM) [14], a pos- 
sibility that has been proposed to solve discrepancies between theoretical predictions and astrophysical observations, 
among which the gamma-ray and microwave emission from the center of our galaxy (eg. lfl5ill6ll . fl7l,[l8il and ref- 
erences therein). It has also been shown that some dark matter interactions could provide an accelerated expansion 
phase of the Universe lfl9l I20I l2lll . (b) The DM could potentially contain more than one particle species, for example 
a mixture of cold and warm or hot dark matter ll22l,l23ll . with or without inter-component interactions. 

In this work we are not concerned with the viability of the different such possibilities, nor with the properties 
of interacting DM models. The aim of this work is to investigate only whether there are repercussions of DM self- 
interactions for the global dynamics of the universe and specifically whether such models can yield an accelerated 
phase of the cosmic expansion, without the need of the dark energy. We note that we do not "design" the fluid 
interactions to produce the desired accelerated cosmic evolution, as in some previous works (eg., j20ll ). but investigate 
the circumstances under which the analytical solution space of the collisional Boltzmann equation, in the expanding 
universe, allows for a late accelerated phase of the universe. 



COLLISIONAL BOLTZMANN EQUATION IN THE EXPANDING UNIVERSE 



As we described in the introduction the DM is usually considered to contain only one type of particles that are stable 
and neutral. In this work, we would like to investigate the cosmological potential of a different scenario, in which the 
dominant "cosmic" fluid does not contain dark energy, is not perfect, and at the same time is not in equilibrium. In 
this context, initially the total energy density is p = Pidm + P» ■ We consider that the self-interacting dark matter does 
not interact significantly with the background radiation, and thus in the matter-dominated epoch, radiation is irrelevant 
to the global dynamics (because of the well-known dependence: p, °< a~ 4 ). Therefore, taking the latter considerations 
into account and assuming that there are no residual radiation products of the DM interactions (otherwise see B24ll ). 
we conclude that in the matter-dominated era the total cosmic dark-matter density reduces to that of the IDM density 
(p ~ Pidm), which obeys the coUisional Boltzmann equation l25ll : 

^ + 3//(f)p + K-p 2 - l P = 0, (1) 
where H(t) = a/a is the Hubble function, *P is the rate of creation of DM particle pairs, and K"(> 0) is given by: 

*=^, (2) 

where a is the cross-section for annihilation, u is the mean particle velocity, and M x is the mass of the DM particle. It 
is obvious that if we impose in the current analysis an effective pressure term such as 

P=(k P 2 -V)/3H, (3) 

then the coUisional Boltzmann equation reduces to the usual fluid equation: p + 3H(p +P) =0. Now if the DM is 
collisionless or the coUisional annihilation and pair creation processes are in equilibrium (ie., P = 0), the corresponding 
solution of the above differential equation is p « a~ 3 , as it should (where a is the scale factor of the universe). In 
contrast, for the case of a non-perfect fluid (ie., having a disequilibrium between the annihilation and particle pair 
creation processes) it becomes clear that we can either have a positive or a negative effective pressure term. 

The possibility of a negative such effective pressure term, ie., the case for which the DM particle creation term is 
larger than the annihilation term (Kp 2 — *P < 0), is of particular interest for its reprecursions on the global dynamics 
of the Universe (see appendix A). In the framework of a Boltzmann formalism, a negative pressure could in general be 
the outcome of dark matter self-interactions, as suggested in IU91 12011 . if an "anti-frictional" force is self-consistently 
exerted on the particles of the cosmic fluid. This possible alternative to dark energy has the caveat of its unknown exact 
microphysics, which is also however the case for all dark energy models. Note, that in the literature, there is a wealth 
of papers discussing the DM particle creation processes in the context of thermodynamics (for example, eg. |26ll ): 

We now proceed in an attempt to analytically solve eq. (1). To this end we change variables from t to a and thus 
eq.ffl]) can be written: 

^=f(a)p 2 + g(a)p+R(a) (4) 

where 

k 3 x V(a\ 

/(«) = Z77 — 7 g(a) = — R(a) = ^ T 1 T. (5) 

aH(a) a aH(a) 

In general it is not an easy task to solve analytically eq. |@), which is a Riccati equation, due to the fact that it is a 
non-linear differential equation. However, eq. (|4|i could be fully solvable if (and only if) a particular solution is known. 
Indeed, we find that for some special cases regarding the functional form of the interactive term, such as *P = ^(a,//), 
we can derive analytical solutions (see appendix B). We have phenomenologically identified two functional forms for 
which we can solve the previous differential equation analytically, only one of which is of interest since it provides 
°< a~ 3 dependence of the scale factor. This is: 

•P(a) = aH(a)R(a) = Q {m + 3)a m H(a) + KC 2 a 2m . (6) 

Although, the above functional form was not motivated by some physical theory, but rather phenomenologically 
because it provides analytical solutions to the Boltzmann equation, its exact form can be justified a posteriori within 
the framework of IDM (see Appendix C). 



The general solution of equation (0) using eq.© for the total energy density is: 



p(a) =Cia m + 
where the kernel function F(a) has the form 

F(a) = exp 



a~ J F(ce) 
[C2-ffx-*f{x)F{x)dx] 1 



-2kC\ 



i H(x) 



-dx 



(7) 



(8) 



We note that kC\ has units of Gyr , while m, C\, and Ci are the corresponding constants of the problem. Obviously, 
Eq.© can be rewritten as 

p(a) = p c {a) + p'(a) , (9) 

where p c — C\ a' n is the density corresponding to the residual matter creation that results from a possible disequilibrium 
between the particle creation and annihilation processes, while p can be viewed as the energy density of the self- 
interacting dark matter particles that are dominated by the annihilation processes. We note that close to the present 
epoch as well as at late enough times (a 1), the p evolves in a similar way to the usual dark matter (see also [27]). 
Finally, if both K and 1* tend to zero, the above cosmological model reduces to the usual Einstein-deSitter model 

(pocfl- 3 ). 

We note that, since p > 0, the constant C 2 obeys the restriction 



C 2 > G(a) = / x~ 3 f(x)F(x)dx > . 



(10) 



Evaluating now Eq.(0 at the present time (a = 1, F(a) = 1), we obtain the present-time total cosmic density, which 
is: p = C] + 1 /C 2 , with d > and C 2 > 0. 



Conditions to have an inflection point 

Using the Hubble function: H 2 (a) = (a/a) 2 = 8nGp(a)/3 , we derive the second derivate of the scale factor: 

& AnGf -2 P - P ^\ . (ID 



a 3 V H 
Inserting Eqs.lfl]) and (01 into Eq.dTTb. we obtain 



a \%G ( Kp 2 -y\ \%G . „ . 

1 P+ u — (p+3P) . (12) 



a 3 V H 

In order to have an inflection point (a = 0), the expression p + 3P = should contain a real root in the interval: 
a £ (0, 1). Thus, with the aid of Eq.©, (01 and ©, we derive the following condition: 

tt " F{ l + .T ° + W^W ~ (m+2)c ' a " H ' ■ <13) 

from which we obtain that m > —2 (where C\ > 0, K > 0, and C2 — G > 0). Evidently, if we parametrize the constant m 
according to m = —3(1 + widm), we obtain the condition widm < —1/3, which implies that the current cosmological 
model can be viewed as a viable quintessence dark-energy look-alike, as far as the global dynamics is concerned. 
We remind the reader that the same restriction holds for the conventional dark energy model in which Pq — wPq 

(w — const < — 1 /3). 
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FIGURE 1. Left panel: The f22,0 ~~ solution space provided by fitting our model to the early-type galaxy Hubble relation of 
Simon et al (2005). Right panel: The corresponding Q.2 o —M x solution space. 



MIMICKING THE A COSMOLOGY 

We will show that for m = the global dynamics, provided by our model (see eqj7j, is equivalent to that of the 
traditional A cosmology. Indeed, we use At = da/(aH) and the basic kernel (Eq.[8]l becomes 



F(a) = exp 



-2kC 



1 



i xH(x) 



dx 



_ e -2KQ (t-t ) 



(14) 



where to is the present age of the universe. In addition, the integral in Eq.(|7]i (see also EqlTOb now takes the form 
G(a) = —KZ(t) and Z(f) = /' a~ i e~ 2KCl ( t ~'°\ We note that at the present time we have G{\) — 0. Obviously, using 
the above formula, the global density evolution (Eq.|7| can be written 



p(a) = C\ + a 



-3 



-2tcCi(t-t Q ) 



\C 2 -G{a)] 



(15) 



As expected, at early enough times (t — » 0) the overall density scales according to p(ct) « a 3 , while close to the 
present epoch the density evolves according to 



p(a)~Ci + - 



a' 



C 2 



(16) 



which is approximately equivalent to the corresponding evolution in the A cosmology in which the term C\ resembles 
the constant-vacuum term (Pa) and the I/C2 term resembles the density of matter (p m ). We note that the effective 
pressure term (Eq.0, for K — > 0, becomes *P ~ ?>C\H, which implies that: P ~ — /3H = —C\. Therefore, this case 
relates to the traditional A cosmology, since C\ corresponds to Pa (see Eq. [Tol l. We now investigate in detail the 
dynamics of the n = model. 

From Eq.(|T5b, using the usual unit-less i2-like parametrization, we derive after some algebra that 



Ho) 



1.0 



^10^2,0" e 



-2KCi(t-f ) 



£2i,o+KCin 2 .oZ(0 



(17) 
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FIGURE 2. Comparison of the Hubble function provided by the traditional ACDM model, which coincides with our m = model 
(for the best-fit model of the two free parameters - see text). The dashed line corresponds to our m = IDM model for the highest 
kC\ bound, provided by our fitting procedure (~ 10~ 3 ). The dot-dashed line corresponds to our K = IDM model for the best-fit 
model parameters (m ~ —0.3 and £22.0 — 0.28). Finally, the points correspond to the observational data ( l28h ). 



where £li q = %kGC\/2>Hq and £^20 = %kG/3HqC2, which in the usual A cosmology relates to £2a and £2 m , respec- 
tively. 

We can now attempt to compare the Hubble function of Eq.(IT7b to that corresponding to the usual A model. To 
this end, we use a % 2 minimization between the different models (our IDM Eq(l7]or the traditional ACDM model) 
and the observed Hubble relation derived directly from early-type galaxies at high redshifts ( ll28lo . For the case of 
our IDM model, we simultaneously fit the two free parameters of the model, i.e., £2^0 ar, d K C\ for a flat background 
(£2i,0 = 1 — ^2,0) with Hq = 72 km/sec/Mpc and frj = H^ 1 ~ 13.6 Gyrs which is roughly the age of the universe of the 
corresponding A cosmology. This procedure yields the best-fit model parameters £22,0 = Jjg and log(jcCi) ~ —9.3 
(with a stringent upper limit ~ —3, but unconstrained towards lower values) where x 2 /d.f. = 1.29 (see left panel of 
Fig. 1). Using Eq.(4) we can now relate the range of values of kC\ to the mass of the DM particle, from which we 
obtain that 

(see also right panel of Fig. 1) and since kC\ is unbound at small values, it is consistent with currently accepted lower 
bounds of M r (~ lOGeV) (see [29] and references therein). The corresponding Hubble relation (Fig. 2), provided by 
the best-fit model free parameters, is indistinguishable from that of the traditional ACDM model, because of the very 
small value of kC\ ~ 10 . For completeness, we also show, as the dashed line, the IDM solution provided by 
M x ~ leV (kC\ ~ 10 -3 ), which is the stringent lower bound found by our analysis. In this case, the predicted Hubble 
expansion deviates significantly from the traditional A model at small a values indicating that it would probably create 
significant alterations to the standard BBN (see lT30tl and references therein). 

Although the present analysis does not provide any important constraints on M x (within our model), we plan on the 
future to use a large amount of cosmologically relevant data to attempt to place stronger M x constraints. 



MIMICKING THE QUINTESSENCE DARK ENERGY 



In this case, we assume that the annihilation term is negligible [k = and f(a)=0 and F(a) = 1] and the particle 
creation term dominates. This situation is mathematically equivalent to the gravitational DM particle creation process 
within the context of non-equilibrium thermodynamics B26I1 . the important cosmological consequences of which were 
studied by [21] (see also references therein). 

Using our nomenclature it is straightforward to obtain the evolution of the total energy density from Eq.©, as: 

p(a) =Da- i + C x a m , (19) 

where D = l/C?. The conditions in which the current model acts as a quintessence cosmology, are given by C\ > 0, 
C2 > 0, and wjDM — — 1 — '«/3, which implies that to have an inflection point, the following should be satisfied: 
widm < — 1/3 or m > —2 (see section 2.1). We note, that the Hubble flow is now given by 

\ =a 2fi a- 3 + a uo a m , (20) 

where Q.20 = 8kGD/3Hq and Q.\ q = 8kGCi/3Hq. Finally, by minimizing the corresponding we find that the 
best-fit model values are Q.20 ~ 0.28 and m ~ —0.30 (widm — —0.90) with ^ 2 /d.f . = 1 .29. The corresponding Hubble 
flow curve is shown in Fig. 2 as the dot-dashed line. We note that this solution is mathematically equivalent to that of 
the gravitational matter creation model of ll2lll . 



CONCLUSIONS 

In this work we investigate analytically the evolution of the global density of the universe in the context of an 
interacting DM scenario by solving analytically the collisional Boltzmann equation in an expanding Universe. The 
possible disequilibrium between the DM particle creation and annihilation processes, regardless of its cause and in 
which the particle creation term dominates, creates an effective source term with negative pressure which (acting as 
dark energy) provides an accelerated expansion phase of the scale factor. Finally, comparing the observed Hubble 
function of a few high-redshift elliptical galaxies with that predicted by our simplest IDM model (m = 0), we find that 
the effective annihilation term is quite small. 



APPENDIX A: THE MAIN ASSUMPTIONS OF OUR MODEL 

The reason why a cosmological constant or a component mimicking it leads to a cosmic acceleration is because it 
introduces in Friedmann's equation, which governs the global dynamics, a component which has an equation of state 
with negative pressure. Our model creates exactly the equivalent to the above, ie., an effective source term which 
has negative pressure. The avenue throw which this is accomplished is via the collisional Boltzmann equation in an 
expanding universe (eg. [25]) in which the disequilibrium between the annihilation (Kp 2 ) and particle creation (¥) 
processes provides the effective pressure: P = (Kp 2 — X i')/3H, which when negative acts exactly as a repulsive force, 
and therefore provides a cosmic acceleration. The basic assumption of this work is that the DM fluid is non-perfect, 
ie., there is a disequilibrium between the particle pair annihilation and creation processes. Note that in our general 
solution (eqJTJl we do not specify the direction of this disequilibrium. Only in the detailed analysis we seek solutions 
which provide an accelerated phase of the universe (for which of course one condition is P < 0), which we indeed find. 



APPENDIX B: SOLUTIONS OF THE RICCATI EQUATION 

With the aid of the differential equation theory we present solutions that are relevant to our eq.©. In general a Riccati 
differential equation is given by 

y =f(x)y 2 +g(x)y + R(x) (21) 

and it is fully solvable only when a particular solution is known. Below we present two case in which analytical 
solutions are possible. 

• Case 1: For the case where: 



R(x) = dnx"- 1 - C\x 2n f{x) - Ci x"g (x) 
the particular solution is x" and thus the corresponding general solution can be written as: 

n -1 



y(x) =Cix n + <f>(x) 



C2 — / f(u)<P(u)du 



where 



<E>(X) = exp 



(2Ciu"f(u) + g(u))du 



(22) 



(23) 



(24) 



andCi,C2 are the integration constants. Using now eq. ©we get y ¥(x) =xH(x)R(x) = C\ (n + 3)x"H(x) + kC 2 x 2 ". 
• Case 2: For the case where: 

R(x)=h'(x) with g(x) = -f(x)h{x) (25) 
the particular solution is h(x) [in our case we have h{x) = —3k~ 1 H(x)]. The general solution now becomes: 



where 



y(x)=h(x)+3>(x) 



<$>(x) = exp 



C2 — / f(u)<P(u)du 



f(u)h{u)du 



(26) 



(27) 



In this framework, using eq. (f5]l we finally get ^(x) =xH(x)R(x) = — 3k 1 xH(x)h' (x). 
Note that the solution of Case 1 is the only one providing a °< a~ 3 dependence of the scale factor (see eqs. [T5l[T6l and 



APPENDIX C: JUSTIFICATION OF THE FUNCTIONAL FORM OF ¥ 

We assume that we have a non-perfect cosmic fluid in a disequilibrium phase with energy density p then from the 
collisional Boltzmann equation, we have that 

y = p + 3Hp + Kp 2 = ^-aH + 3Hp + Kp 2 . (28) 
da 

Furthermore, we assume that for a convenient period of time, the cosmic fluid, in an expanding Universe, is slowly 
diluted according to p ~ C\(X m (m < 0). From a mathematical point of view, the latter assumption simply means that 
a solution of the form °< a" 1 is a particular solution of the Boltzmann equation. Therefore, we have finally that: 

¥ ~Ci(m + 3)a'"H+KC 2 a 2m . (29) 
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